Theoretical analysis of quantum dynamics in 1D lattices: Wannier-Stark description 
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This papers presents a formalism describing the dynamics of a quantum particle in a one-dimensional 
tilted time-dependent lattice. The description uses the Wannier-Stark states, which are localized in 
each site of the lattice and provides a simple framework leading to fully-analytical developments. 
Particular attention is devoted to the case of a time-dependent potential, which results in a rich 
variety of quantum coherent dynamics is found. 

Pacs number(s): 03.65.-w, 03.75.-b, 32.80.Lg, 32.80.Pj 



I. INTRODUCTION 

Quantum dynamics in a periòdic lattice is one of the 
oldest problems of quantum mechanics, whose basis have 
been settled by Bloch and Zener wM, in the 30 's. Aimed 
at the description of the electron motion in crystalline lat- 
tices, this problem has largely been considered, for about 
half a century, as an acadèmic one, because dissipation 
cffccts forbid the observation of most quantum effects 
in the motion of a crystalline electron. Làser cooling of 
àtoms has brought a revival of the interest on such prob- 
lems, as it produces àtoms whose de Broglie wavelength 
is comparable to the wavelength of the light interacting 
with the àtoms, and whose resulting kinetic energy is 
comparable to the typical lightshift induced by the radi- 
ation. The latter feature means that the cold àtoms can 
be trapped in light potentials (or dipole potentials). The 
former means that the atom dynamics in such a potential 
is, in absence of dissipation, essentially quantum. More- 
over, the main source of dissipation is spontaneous emis- 
sion, which can be arbitrarily reduced (if one disposes of a 
powerful enough làser), whereas keeping a constant light 
otential, just by an increase of the laser-atom detuning 



Light potentials are a consequence of the displacement 
of atòmic levels resulting from the interaction with light, 
corresponding to a process in which a photon is absorbed 
transferring the atom to an (virtual) excited state from 
which the àtoms de-excites back to the ground-state by 
stimulated emission. Such a process induces an energy 
shift of the atòmic levels that can be deduced from sec- 
ond order perturbation theory and which is proportional 
to light intensity, to the square of the coupling (that is, 
to |d eg • e| 2 , where d eg is the dipole matrix element be- 
tween the states g and e, and e is the polarization vector 
of the light), and to the inverse of the laser-atom detun- 
ing Sr, = — uj eg (uj eg is the Bohr frequency between 
states g and e). Any spatial gradient of this energy shift 
produces a (conservative) force, and thus a potential. A 
simple example is that of a standing wave formed by two 
counter-propagating parallel-polarized beams. Placed in 
such a standing wave, an atom perceives a periòdic one- 



dimensional potential whose strenght varies sinusoidally 
in the space. Standing waves (with little variations) form 
the model potential considered in the present work. 

Light potentials generated by standing waves have 
been used in many experimental studies of quantum 
dynamics. For example, Bloch oscillations have been 
observed both with single àtoms Q and with a Bose- 
Einstein condensate (BEC) || in an accelerated standing 
wave. Wannier-Stark ladders j(| and collective tunneling 
effects JT) have also been studied with such a system. 
Àtoms placed in an intense, phase-modulated, or pulscd, 
standing wave realize a paradigmàtic system for theo- 
retical and experimental studies of quantum chaos, the 
so-called Quantum Kicked Rotor fl8|-{l·l|] . 

In this paper, we consider the quantum dynamics of 
an atom (of mass M) placed in a tilted sinusoidal poten- 
tial whose phase (that is, the position of its nodes) can 
be modulated in an arbitrary way, corresponding to the 
Hamiltonian 



(i) 



where Xo(t) is a phase and /(£) a force, both being 
(eventually) time-dependent, and = 2t:/\l is the 
wavenumber of the standing wave. Different temporal 
dependences of xo(t) can be considered. For instance, the 
accelerated case, Xo(t) = (l/2)aí 2 , which has been stud- 
ied in 0^,0], is equivalent to an inertial force F = Ma 
in the frame of the potential (see appendix |Aj) . 

A natural energy unit in such a context is the "recoil 
energy" , defined as the change in kinetic energy of the 
atom corresponding to the absorption a photon, given 
by 



2 1,2 



Er 



2M 



(2) 



to which one can associate a recoil frequency lür = Er/U 
and a recoil momentum pr = hkr,, etc. It is also useful to 
re-scale the variables: X = x/(Xl/2), where Al/2 is the 
step of the periòdic lattice, and r = lúrí. With these def- 
initions the above Hamiltonian takes the following form, 
which is retained in the rest of the paper: 
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p2 

H= — + V cos{2Tr[X-X (T)}}+F(T)X (3) 

where Vo = v /E r , F = J\l/2Er (note that in this 
system the momentum operator in the real space is P = 
—í(d/dX), h = 1, the reduced mass is m* — tt 2 /2 and 
d = 1 is the step of the lattice). For simplicity, in what 
follows, we shall write the rescaled variables x, p and t. 

In the next section, we briefly review the time- 
independent Hamiltonian case in Eq. (^), introduce the 
Wannier-Stark basis, and show that it leads to a very 
simple description of the Bloch oscillation. In the fol- 
lowing sections, we shall discuss the more complicated 
dynamics that arises when a harmònic modulation of the 
lattice is applied. 

II. THE BLOCH OSCILLATION IN THE 
WANNIER-STARK DESCRIPTION 



resonances) have long enough life-times compared to the 
experimental times. The existence of these states has 
been evidenced in 1988 in a semi-conductor superlattice 
|p^| , and 1996 with cold àtoms in an optical lattice M. 
Note that, being spatially localized, Wannier-Stark states 
provide an ideal tool for the description of the wave func- 
tion of a cold atom (as produced by a "Sisyphus-boosted" 
MOT) placed in the potential, whose de Broglie wave- 
length (around Ai/3) is of the order of the lattice period 
Al/2. 

Consider the properties of the time-independent 
Hamiltonian Hq 

H a = -ï— + V Q cos(2ttx) + Fx (4) 
2m* 

where F is a constant force. The eigenfunctions of H 
are Wannier-Stark states forming an energy ladder whose 
separation is = Fd, E nm = E m + nFd. 



The Bloch oscillation (BO) is a well known phe- 
nomenon discovered by Zener while studying the quan- 
tum properties of an electron in a (perfect) crystal sub- 
mitted to a constant elèctric field The BO arises 
when a small spatial tilt is added to the lattice: Quan- 
tum partides do not fall along the slope of the potential, 
but perform a periòdic, space-limited, oscillation. BO is 
thus a strictly quantum behavior. We shall call "lattice" 
the untilted potential, and "tilted lattice" the sum of the 
lattice and the tilted potential. 

The BO is usually described in the basis of the so- 
called Bloch states i.e., the eigenstates of the lattice. 
In this paper, we use another framework corresponding to 
the eigenstates of the tilted potential. While the lattice 
is invariant under spatial translations by a múltiple of 
the spatial period d, the tilted lattice has a more compli- 
cate symmetry: it is invariant under simultaneous spatial 
translation by a lattice period d and energy translation 
by lüb = Fd (ujb is the "Bloch frequency" ). One then ex- 
pects the eigenenergies to form "ladder" structures sep- 
arated by wb, the so-called Wannier-Stark ladders, in- 
troduced by Wannier in connection with the problem of 
electrons in a crystal submitted to a homogeneous elèc- 
tric field |^] . Each element of the ladder corresponds to 
eigenfunctions (the Wannier-Stark states) centered at a 
given well, and thus separated by an integer múltiple of d. 
The form, and even the existence of these eigenstates has 
been the object of a long controversy, that has been set- 
tled only recently |Ï3]| . In the present paper, we consider 
a spatially limited lattice, extending over many periods, 
and limited by an infmite-height box. This changes only 
very slightly the "bulk" properties of the system, and the 
eigenenergies and eigenstates obtained numerically dis- 
play (to a very good approximation) the expected ladder 
structure described above. In this framework, Wannier- 
Stark states (WSS) are the eigenfunctions of the system. 
The relation to the case of an infinite tilted lattice makes 
no problem if the corresponding states (Wannier-Stark 
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FIG. 1. (a) Periòdic potential with a tilt. (b) The WSS 
state íp n is localized in the well of index n = 1 and has ap- 
preciable overlap with neighbor lattice sites n = 0, 2. This 
eigenfunction is obtained for Vo = 2.5 and F — 0.5. The 
"numerical" box includes 64 lattice sites. 



The BO can be advantageously described by using the 
Wannier-Stark state (WSS) localized inside a given in- 
dividual lattice well, corresponding to the lowest energy 
of the states associated to this well (see Fig. |^). Note 
that considering only the ground state of each well is 
equivalent to the restriction to the first Bloch band in 
the description based on Bloch states. We also choose 
strong enough F and Vq to produce well-localized WSS 
M|. The WSS associated to the lattice well labeled n 
is noted tp n {x) (supposed real) and the corresponding 
eigenenergy is E n (we drop the index m). The symmc- 
tries of the potential discussed above then imply: 

(p n+P (x) = <p n {x - pd) (5) 

and 

E n+p = E n +pFd. (6) 

Let us describe the atòmic wave function by a super- 
position of WSS 



2 



V(x,t) =^2c n (t)(p n (x) (7) 

with c n (t) = c„e~ ***; c n being the amplitudc at t = 0. 
The dynamical quantities of the atom can easily calcu- 
lated. For instance, the mean value of the atòmic position 
operator is: 

(X) =J2 X ^n K\ 2 + ]T ( X n,mC* n Cme í(En - Em)t + C.c) 
n n <m 

(8) 

where X niTn = (ip n \x\ip m ). As long as the WSS is 
well localized in the respective well, we can keep only 
nearest-neighbors contributions (n = m ± 1) and derive 
a simplified expression: 

(x) = 5 + X„ >n+1 ^ c* n c n+1 e-^ Bt + c.c^J (9) 

where x = ^2X n ^ n \c n \ 2 is the mean position and 
X nn+ i = Xq^i = Xo,-i = X n .„_i is independent of n 
fÏ6f . This result evidences a counter-intuitive property of 
the quantum motion in a tilted lattice: instead of falling 
along the slope, the atom performs an oscillation with 
frequency lüb (the Bloch frequency). The amplitude of 
this Bloch oscillation is proportional to Xo.i and grows 
with the overlap between neighbors WSS, i.e for a small 
slope F and small lattice depth Vq. The physical origin 
of the BO appears here clearly as an interference effect 
between neighbor sites, as c*c„+i is the coherence be- 
tween the sites n and n + 1. Note the lack of oscillations 
if the atom is localized in only one well. Eq. (||) also pre- 
dicts the occurrence of harmonies with frequency ptos (p 
integer) but with smaller amplitudes since they involve 
the coupling strength X nt7l+P [|Ï7f . 

The description of the BO given here is very different 
of the usual "solid-state" approach. There, the Bloch 
states (eigenstates of the untilted lattice) are taken as the 
reference basis. The oscillation is described in a semi- 
classical frame as the periòdic evolution of the atom's 
quasi-momentum Ft (in the first Brillouin zone) with 
period T = 2n/(Fd) |Ï|Q. Although intuitive, this 
approach does not make clear the rolc of quantum inter- 
ference produced by the periòdic lattice structure as the 
bàsic mechanism underlying the Bloch oscillation, which 
is evidenced in our approach. 

III. THE MODULATED POTENTIAL IN THE 
WANNIER-STARK DESCRIPTION: RESONANT 
DYNAMICS 

With the existence of the natural frequency ujb of the 
system in mind, one is tempted to investigate the quan- 
tum dynamics in presence of a harmònic external forcing. 
The WSS approach proves to be very efficient, since it 



allows a fully-analytical description. After some general 
considerations, we study in this section the case of res- 
onant forcing, and show that it leads to a very rich and 
interesting dynamics. The general (non-resonant) case 
will be treated in the next section. 

Consider the Hamiltonian of Eq. (||) with a constant 
force F and a lattice phase modulation 

xo(t) = asm(íjí). (10) 

The developments are simpler if we use a unitary trans- 
formation that transforms the modulation in a time- 
dependent force. Physically, this is equivalent to move 
to an accelerated reference frame in which the lattice is 
at rest adding thus an inertial force (AppendixH). 

In this frame, the new Hamiltonian is given by Eq. (0) 
plus a time-dependent force F'(t): 

d 2 x (t) 

F'(t) = m* — - — = — m*aLü 2 sin(wí) = — Fosin(wí) 

(11) 

where Fq = m*aui 2 is the amplitude of the inertial force. 
The harmònic time-dependence in Eq. (O) is the analog 
of an AC elèctric field for electrons in a (perfect) crystal. 
The dynamics in such a system can be described in a 
quite simple fashion by writing the state of the atom as 
a superposition of WSS, Eq. ([?]). The coeficients c„(t) 
can be obtained by reporting Eq. (Q) into the Schròdinger 
equation 

[H - F x S m(ujt)}*(x,t) = i^A^ 1 (12) 

at 

where H is given by Eq. (^), with eigenstates ip n (x). 
This produces the following set of coupled differential 
equations for the c n (t): 

c n (t) = -iE n c n {t) + iF sin(wí) ^2 X n (t) (13) 

in 

where c n = dc n /dt. Neglecting temporarily the coupling 
between different WSS, (i.e putting X„ jTn = J n 
the amplitudes are obtained simply as c„ = cxp[iç!> n (í)], 
with the time-dependent phase 

4> n (t) = -E n t - cos(wí) (14) 

Lü 

where X n;n — Ao.o + nd, depends on the site index n 
@|. We now write c n (t) = d„(í)e^"W. The amplitudes 
d„ obey the following system of differential equations: 

d n = íFq ^ X n ^ m d m (t) exp{i [<p m (t) - 4> n (t)}} sm(ut) . 

(15) 

After Eq. (|Ï4|), the phase difference <fi m (t) — <fi n (t) is: 
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F d , .' 
lobí H cos(o;í) 

Lü 



4> m {t) - 4>n{t) = (n - m) 
where we used Eq. (|). Eq. |TJ) can be recast as 



(16) 



d„ = iF ^2 X P d 

n+p 



sm(cjí) 



(17a) 



{ 



p^O i x 

3 i[(( + l)tJ— pUB]t P Í\{1— pusjí 



} (17b) 



where, X p = X„. n+p [JÏ6j. J n (x) is the Bessel function of 
the first kind, and we have used the well-known property 
of Bcsscl functions: 



+ OC 



-izcos(ut) 



Uíüt 



(18) 



Z= — oc 



Note that the sum over lattice sites (i.e over p) extends 
only over a few neighbors sites, since the coupling coefi- 
cients X p rapidly shrink to zero [|Ï7f . If the modulation is 
smooth enough to avoid projections on other states of the 
system, the sum over the harmonies of the modulation 
(i.e. I) is also limited to a few terms closc to / = (typ- 
ically, l ni ax ~ pFod/uj = pm*aujd ~ 0(1))- Therefore, 
only a finite number of terms are to be retained in the 
above expression. On th e other hand, the evolution of d n 
described by Eqs. ( 17b ) is a sum of oscillations with fre- 
qüències [{l ± Í)lu — pus]- In the following we keep only 
the so-called secular terms, that is, terms that oscillate 
slowly or do not oscillate at all. The resulting "close to 
resonance" dynamics is observed when (l ± l)u> w plüb, 
i.e when the forcing frequeney lo is commensurable (or 
almost) with the system's natural frequeney lüb- 

Let us consider the simpler resonant case, u> = tog. 
Due to the relative strength of the factors X p we can 
keep, to a good aceuracy, only the contribution of the 
next-neighbor site (p — 1), which leads to the following 
expression: 



d n (t) = fli [d n+í - dn-x] 



(19) 



where 



2 [ \ Lüb J \U B 

üj b X 1 ÍF d\ 

This equation is similar to a "dipole coupling" between 
sites n and n ± 1 where Í2i plays the role of a Rabi fre- 
queney. Note that, contrary to intuition, the coupling 
towards the left or towards the right neighbor is the same. 

The meaning of Eq. ( |Ï9| ) can be better appreciated by 
searching for the plane-wave solutions of the form: 



dn(t) 



^i(kç,dn-\-Lüt) 



(21) 



The (dimcnsionless) wavenumber ka takes into account 
the phasc diffcrence between neighbor sites. Substitution 
into Eq. (|Ï9| ) leads to the dispersion relation 



lü = 2í7i sin(fcod) 
and to the group velocity v g = dej/ dko 
v g = 2ílidcos(kod) . 



(22) 



(23) 



This result shows that the dynamics is different depend- 
ing on the wave number fco. For instance, if k^d = ±7r/2 
(phase-quadrature from site to site), v g — and there 
is no global motion. If kod — tt, v g = —2Q.id and the 
global motion is a fall along the slope of the potential with 
the maximum speed 2íl\d. More interesting is the case 
ka = 0, where v g = 2Çl\d: the atom then climbs up the 
slope of the potential with a constant maximum speed: 
there is, in this case, coherent transfer of energy from the 
modulation to atom, thanks to the particular phase re- 
lations between neighbor sites. Note also that, contrary 
to the motion of a classical particle, the speed \v g \ is in- 
dependent of the sense of displacement: the wavepacket 
climbs the slope up or down at the same speed. 

More detailed information on the wavepacket motion 
can be grabbed by writing the amplitude d n in the more 
general form: 



d n(t) = f n (t)e 



i(kodn-\-ujt) 



(24) 



where f n are complex amplitudes describing the enve- 
lope of the atòmic wavepacket, assumed to vary slowly 
in time as compared to the frequeney ljb, an d in space 
as compared to the lattice period d. Reporting the above 
expression into Eq. (pa), we get: 



/„ + iwf n = fil [C0s(k a d) (fn+1 - fn-x) 

+*sm(*cod)(/»+i + /n-i)] ■ (25) 



Using the dispersion relation Eq. (|22|) and keeping only 
slowly varying contributions, we obtain: 



/„ = fli [cos(k Q d) (fn+i - fn-x) 

+ÍSÍXí(k d) (fn+X +fn-X~ 2/„)] . 



(26) 



Since f n vary slowly in space, we can take the continuous 
limit (with respect to the variable x — nd) and deduce 
the following equation for the wavepacket envelope: 

f{x,t) = fli ^dcos(k d)—+td 2 sin{k d)-^)f(x 1 t) . 

(27) 

This equation is an interesting piece of information. If 
kod = ±7r/2, one has a diffraction equation of the form 
/ = ±iflid 2 d^f , which describes the spreading of the 
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atòmic wavepacket (i.e. diffraction) with no global dis- 
placement. The case kgd = 0,n gives the wave equation 
/ = Vgd x f (v g — 2dQi), describing a wavepacket trav- 
cling with constant velocity v g and no deformation: the 
wavepacket presents an ascending or descending coherent 
motion. Mixed behaviors, i.e. spreading at a diffusion 
rate Q,\d 2 sin(fcod) and uniform displacement with group 
velocity v g = 2Çl\d cos(fcod), are found for other vàlues 
of fco- The general solution can be easily obtained by 
performing a spatial Fourier transform of Eq. (|27j) . 
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FIG. 2. The atòmic wavepacket is obtained from a full in- 
tegration of the Schrodinger equation corresponding to the 
modulated tilted lattice. The wavepacket at t = has a 
site-to-site phase of n/2. The wavepacket spreads in time 
with no displacement, in an agreement with the theoretical 
prediction. Parameters are Vo = 2.5, F — 0.5, u — lüb = 0.5 
and a = 0.2 
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FIG. 3. Same as Fig. ^, except that the site-to-site phase 
is ko — 0. The wavepacket has a shape-preserving motion up- 
wards the slope of the potential. The observed group velocity 
is 0.030, in good agreement with the prediction of Eq. (p3|), 
v g = 0.032 (X ,i = 0.13). 



In Fig. ||, the wavepacket, prepared with ko = has 
a uniform displacement while preserving its shape. Its 
group velocity obtained from the numcrical simulations 
is Vg = 0.030, in very good agreement with the theoretical 
value from Eq. (03) which is v g — 0.032. 
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FIG. 4. Evolution of a wave packet for ui = 2lüb- The ini- 
tial wavepacket is a superposition of ip n (n odd) which are in 
phase, and of ip„ (n even) with phase difference ir. The initial 
wavepacket separates in two packets with opposite velocities. 



Other rcsonant behaviors are observed if uj = quB (q 
integer). From the general expression of Eq. (|Ï7|), onc 
finds for instance a next-to-neighbor (n — > n±2) resonant 
interaction if lo — 2lüb, leading to: 



d n (t) = íl 2 [d n +2 - d n - 2 ] 



(28) 



with íl 2 = (cü B /d)X 2 Ji(F d/L·JB)· Note that X 2 < X 1 
JÏ7t . We illustrate in Fig. [| the temporal evolution of 
a wavepacket, obtained by numerical integration of the 
Schrodinger equation. The initial state is prepared as a 
superposition of two packets: one packet is constructed 
with in-phase amplitudes (that is, c n and c„+2 have the 
same phase for n odd), and the other one is constructed 
with amplitudes in phase-opposition (that is, c„ and 
c n+ 2 have opposite sign for n even). The first packet 
moves with velocity v g — 4ÍÍ2CÍ, and the second with 
v g = — The figure displays an original behavior 

showing each of these two initially inter-penetrated pack- 
ets moving independently in opposite directions, creating 
a highly delocalized state. 



IV. THE MODULATED POTENTIAL: GENERAL 
CASE 



The above result shows that, depending on the initial 
wavenumber kg, i.e. on the way the initial wave packet 
is prepared, the atom behaves in qualitatively different 
ways. Figs. ^| and ^ are obtained by a direct integration 
of the Schrodinger equation, with the Hamiltonian given 
by Eq. (|ï^), for different evolution times. In the first case 
the initial wavepacket is prepared with site-to-site phase 
kgd = 7r/2 and diffraction is clearly visible, as predicted. 



In this section we generalize the results of the preceding 
section to the case of a non-resonant modulation. We 



follow essentially the same steps as in Sec. III, and we 
shall skip al gebr aic details of the calculations. Corning 
back to Eq. (17a) and looking for a solution of the form: 



d n {t) 



= lp *! + </>(*)] 



(29) 



one gets the instantaneous frequency: 



5 



4> = 2F ^2 x p cos{p[k Q d - 6(t)]} sin(uí) (30) 

p>0 

where 0(t) = Wst + (F d/uj) cos(wí), and we used X p — 
X-p. The group velocity v g = d(f>/dko is thus 

v g = 2F d^2pX p sm{p[6(t) - k d}} sin(wf) (31) 

As in the preceding section, more detailed behavior is 
obtained by putting 

d n (t) = f n (t)e^' dn +^ (32) 

where /„ are slowly-varying amplitudes. The generaliza- 
tion of Eq. ( p6| ) is then: 

/„ = iF Q XAfn+v - fn}e tp{kod - e) sin(wt) (33) 

or, 

fn = Fo^Xp sin(cjí) 

p>0 

{- [fn+p - fn-p] sm[p(k a d - 9)} 

[(fn+p + fn-p - 2/ n )] CO S [p(k d - 6)}} . (34) 

Taking the continuous limit of the above expression then 
produces an equation describing both the propagation 
and the diffraction of the wavepacket: 

(d d 2 \ 

f(x, t) = [v g (t)— + mt)g^J f(x, t) 

+2iF X p cos[p(é> - k d)] sm(ut)f(x, t) (35) 

p>2 

where 

D{t) = F d 2 J2p 2 X p cos[p(6 - k d)\ sin(wí) (36) 

p>0 

and the group velocity v g is given by Eq. @. Note 
that the last term in Eq. (|35|) is a phase term which is 
0(X2) « 1 and does not contribute to the probabil- 
ity density \f(x,t)\ 2 . For the sake of lightncss, it is not 
considered in the following. 

The Fourier transform of Eq. (|35| ) with respect to x 
produces an algebraic equation for the Fourier transform 
f(k,t) of f(x,t), whose solution is: 

/(fe, t) = e ife 'We ifc3A W/(fe, 0) (37) 

and thus: 

f( X) t) = -L· f e ik[x+x'(t)] e ik 2 A(t) ( 38 ) 
V 2n J 

with 

x'(t) = / Vg(r)dT (39) 
Jo 



and 

A(i) = f D(r)dT . (40) 
Jo 

This expression describes a coherent motion of the 
wavepacket formed of an oscillatory motion with the 
time-dependent group velocity Eq. (|3l"|), and a diffusive 
motion with a time-dependent diffusion coefficient D(t). 
For example, if one builds an initial gaussian packet of 
width ao, f(x,0) = exp(— x 2 /ag) one finds, after some 
straightforward calculations: 



where 




Position x 

FIG. 5. Spatio-temporal behavior of ^(a;, í)| 2 obtained nu- 
merically by integration of the Schrodinger equation (gray 
level convention: the maximum vàlues of |Ç(a;,í)| 2 are de- 
picted in black), from í = to í = 4n / 'S (that is, 2 periods 
of the beat frequency S = u) — ujb = 0.02). Other parameters 
are the same as in Fig. 



The physical meaning of our developments can be ev- 
idenced by considering the case where cu differs from lüb 
by a small detuning S — lü — lüb,\S\<^ujb, and keep only 
leading-order terms of order 0((5 _1 ). The wavepacket 
then undergoes a harmònic oscillation at the beat fre- 
quency S with a group velocity given by Eq. ( |3Ï] ) 

v g {t) = 2ü 1 dcos(k d + St) (43) 

corresponding to a periòdic mean position displacement 



G 



(x(í)) = x(0) - [sin(fc d + St) - sin(fc c/)] (44) 
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íl d 2 

A(t) = — - — [cos(k d + St) — cos(k a d)] 



(45) 



The results of Sec. [II for a resonant excitation are nat- 
urally recovered in the limit 5 — > 0. 
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APPENDIX A: UNITARY TRANSFORMATION 



FIG. 6. Mean vàlues (a) (x(t)) and (b) y/(x 2 (t)) - (x(t 
obtained from the wavepacket dynamics depicted in Fig 



We have performed the integration of the Schròdinger 
equation for a detuning S = 0.02. Fig. ^| shows the 
spatio-temporal dynamics of the wavepacket and clearly 
evidences the periòdic oscillations and the breathing at 
frequency S predicted above. Figure [| shows the mean 
position and width of the wavepacket as a function of 
time. The comparison with the theory is very good. For 
instance, the amplitude of (x(t)) is found numerically as 
1.55 compared to the theoretical value 1.73. 



Eq. ( |11| ) is obtained if we perform a unitary transfor- 
mation 



U(t) 



= e ix (t)p e -iP(t)x e *7(*) 



(Al) 



where we have included translation operators in space 
and momentum with (3 = 2m*Ao(i.e momentum of a 
particle of mass m*). In this framework, following p0[ , 
we obtain (with UXU+ = X + X and UPU+ =P + j3) 



2m* 



H' = U HU + + i (d t U) U+ 
+ V cos(2ttA) + F(X + X (t)) - X (t)P + 

${t) (A + A )-7 



with 7 = FX + m* Xo A + (m*/2)Xg : 



H' = 



Jz 

2m' 



+ Vo cos(2ttA) + {F + m* X (t))X (A2) 



Therefore, in the frame of the periòdic potential, the 
Hamiltonian contains an inertial force proportional to 

Xo (t). 



V. CONCLUSION 



We have studied, in a fully-analytical way the dynam- 
ics of a wavepacket in a static and time-modulated tilted 
potential in the framework of the Wannier-Stark states. 
This basis is well suited to the description of the state of 
a cold atom (as produced by a Sisyphus-boosted MOT). 
Moreover, it provides a simple description the atòmic dy- 
namics, which is proved to be very rich: A variety of co- 
herent motions are obtained depending on the prepara- 
tion of the initial wavepacket and its site-to-site quantum 
coherence. We can note that the description introduced 
here is, in its principle, independent of the details of the 
lattice, provided it presents localized states in the lat- 
tice sites. It is therefore generalizable to other kinds of 
lattices. 

Let us finally, mention that the present work distin- 
guishes from the more usual "solid-state" approach which 
is based on Bloch functions. We postpone for a fortheom- 
ing work the detailed comparison between the Wannier- 
Stark and Bloch approaches. 
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